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ABSTRACT 

Aims. We study the deuteration and ionization structure of the DM Tau disk via interferometric observations and modeling of the key 
molecular ions, HCO^ and DCO^. 

Methods. The Plateau de Bure Array is used to observe DM Tau in lines of HCO^ (1-0), (3-2) and DCO+ (3-2) with a ~ 1.5“ angular 
and ~ 0.2 km s“* spectral resolution. Using a power-law fitting approach the observed column densities profiles are derived and thus 
the isotopic ratio = DCO^/HCO^. Chemical modeling allowed an exploration of the sensitivity of HCO^ and DCO^ abundances 
to physical parameters out with temperature. A steady state approximation was employed to observationally constrain the ionization 
fraction x(e^). 

Results. Fitting of radiative transfer models suggests that there is a chemical hole in HCO^ and DCO^, extending up to 50 AU from 
the star. More work is required to discern the cause of this. The observed column densities of HCO^ and DCO^ at 100 AU were 
(9.8 (;q 7 )x 10'^ and (L2±0.7)x lO'^ cm“^ respectively. Where both HCO^ andDCO^ were present, Ro was found to increase radially 
from 0.1 at 50 AU to 0.2 at 450 AU. This behaviour was well reproduced by the chemical model. The X-ray luminosity of the central 
star, the interstellar UV and CO depletion were found to be the most important physical parameters controlling the abundances of 
HCO^ and DCO^. Differences in the vertical extent of HCO^ and DCO^ molecular layers resulted in different responses to changing 
physical parameters, manifesting as radial gradients in R^. The ionization fraction was found to be x(e^) ~ 10“^ in the molecular 
layer, comparable to the disk averaged value. Modeling shows that while HCO^ is the most dominant charged molecular ion in our 
disk model, atomic ions, such as C^, S^, H^, Na^ and Mg^, dominate the charge in both the molecular layer and disk atmosphere. 
Conclusions. A high value of Rd is indicative of continued deuterium fractionation in a protoplanetary disk after pre/protostellar 
phases. Radial properties of can be employed to discern the importance of ionization from X-rays and UV, thus necessitating 
the need for more, high resolution observations of DCO^ and other deuterated species in disks. A steady-state approach commonly 
adopted for constraining ionization degree in prestellar cores is not applicable for disks where accurate determination of the ionization 
fraction in the molecular layer requires knowledge of the atomic ions present as molecular ions are relatively sparse. 
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1. Introduction 


In the view of recent exciting discoveries of various extraso¬ 
lar planets with the Kepler satellite, including Earth-like planets 
in a habitable zone ([Quintana et al.|2014|l, and recent ground¬ 


breaking Subaru, Hers chel, and AL 


toplanetary disks (e.g., jOberg et al. 


Williams & Cieza|20 11 [ Eedele et al. 


VIA observations of pro- 
2010 jSturm et al.||2013 


2013a 


Grady et al. 2013 


van der Marel et al.||2013[ [Rosenfeld et al.||2014] l, it is with 
great anticipation that we begin to unravel the planet formation 
process. Interferometric observations of protoplanetary disks 
(PPDs) provide a crucial tool in the quest to understand these 
enigmatic objects. 

The high-resolution measurements of the dust continuum 
and molecular emission lines of various optical thicknesses al- 


* Based on observations carried out with the IRAM Plateau de Bure 
Interferometer. IRAM is supported by INSU/CNRS (France), MPG 
(Germany) and IGN (Spain). 


low for probes of physics and chemistry in distinct disk regions 
(for recent reviews, see Henning & Semenov 2013[ [Quintana' 
[et al.|2014) l. To make sense of these data one has to employ disk 
physical/thermo-chemical models that predict the thermal and 
density structures of disks (e.g., [Woitke et al.[[2009l [Akimkin 
[eraL][20T3| ), gas-grain chemical models to calculate molecular 
distributions (e.g., Willacy et al.|1998 Aikawa & Herbstj 1999b 


[Semenov et al.[[2010[ [Walsh et al.[[2010[ ), and radiative transfer 


3Ep 

Hogerheijde & van der Tak 

2000 

van Zadelhoff et al. 

2002 

Pinte et al. 2006[ Pavlyuchenkov et al. 2007[ Jprgensen 

et al. 

2014|l. Another approach is to perform an iterative fitting of 


the observed spectra/interferometric visibilities, applying sim¬ 
pler power-law models of disk physical and molecular structure 


* http://youngstars.nbi.dk/artist/Welcome.html 
^ http: //WWW. ita.uni-heidelberg.de/~dullemond/ 
software/radmc-3d/ 
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coupled to a fast LTE/LVG radiative transfer model (e.g., Guil 


loteau & Dutrey|1998| Pietu et al.|2007[ Qi et al.|2008[ Rosen- 

feld et al.|2013|l. A combin 

ation of these approaches have been 
our “Chemistry In Disks” (CID) con- 

used in series of articles by 

sortium (Dutrey et al.|2007[|Schreyer et al.|2008 Henning et al. 

2010 Semenov et al.|2010 

[Dutrey et al.|201 1[ Guilloteau et al. 


(|Harju et aL||2006 

Hogerheijde et al.| 

20061 Parise et al.||2011t 

[Vastel et al.|[2006 

Parise et al. 

2011 

1 and have not yet been 


|lon et al.|201 1 1 . Thus, to fully characterize deuterium chemistry 


Table 1: Stellar and disk properties of DM Tau with values taken 
from Dutrey et al. (2007 1 and Henning et al. (2010 1 . 


2012 I. 

In this paper we present interferometric observations of 
HCO+ J = (3-2), J = (1-0) and DCO+ J = (3-2) in DM Tau, 
a bona fide T-Tauri star in the Taurus-Auriga star-forming re¬ 
gion. Owing to its large size, ~ 800 AU, and moderate inclina¬ 
tion, i 35°, it has become one of the best studied PPDs. In 
one of the first millimeter studies of PPDs, [Guilloteau & Dutrey 
( 1994| l used the '^CO and '^CO J = (2-1) lines to derive a disk 
mass of 1.4 X IO^^Mq and disk radius of =» 700 au. Later, high 
spectral- and spatial-resolution observations of ^^CO, '^CO and 
C'^O have allowed the detection of a vertical temperature gradi¬ 
ent within the disk and the presence of very cold, ~ 10 - 15 K 
CO gas ( Dartois et al.||2003] l, a better constrained disk mass of 
0.05 Mq and an outer radius in CO of 800 AU. 

|Pietu et ^ ( |2007| l used higher resolution observations of 
these lines supplemented with HCO^ J = (1-0) to better constrain 
the disk physical structure and kinematics in the radial direction. 
They confirmed the presence of the vertical gas temperature gra¬ 
dient in the DM Tau system, found that CO has an extended 
distribution in vertical direction, and that the slope of the CO 
surface distribution changes its value with radius. 

In the first CID paper of [Dutrey et al.| P007| l, a sensitive 
observation of N 2 H^ and HCO^ towards three disks (includ¬ 
ing DM Tau) was performed with PdBI, followed by advanced 
physico-chemical modeling. It was found that HCO^ is a ma¬ 
jor polyatomic ion in disks, and that its column density agrees 
with the modeled values at an evolutionary stage of a few mil¬ 
lion years. The ionization degree in the HCO^ molecular layer 
was also derived, ~ 2 x 10“®. 

In the next CID paper by jSchreyer et al.| ( |2008| l the chemical 
content of the DM Tau disk was compared to the disk around a 
hotter Herbig AO star; AB Aur. We found that while the AB Aur 
disk possesses more CO, it is less abundant in other, more com¬ 
plex molecular species compared to the DM Tau disk. This find¬ 
ing gives a hint that high-energy radiation from the central star 
may be important not only for disk thermal structure but also 
for its chemical com plexity (see also Fedele et al.|201 1 2013b 
Oberg et al.j201 la i. 


DM Tau Stellar and Disk Properties 

Right Ascension (J2000) 

04*’33“48';733 

Declination (J2000) 

-h18‘’10'09'.'89 

Spectral Type 

Ml 

Effective Temperature (K) 

3720 

Stellar Luminosity (Lq) 

0.25 

Accretion Rate (Mo yr“') 

2 x 10^’ 

Disk Mass (Mq) 

0.05 

Rout (au) 

800 


Aikawa & Herbst|1999al|Albertsson et al.|2013l|Ceccarelli et al.| 

2014) 1. 


The first detection of a deu terated species in a disk was made 

who detected DCO^ in TW Hya 


by van Dishoeck et al. (2003 


and discerned a disk averaged ratio of Rq ^ 0.04, a value similar 
to that found in pre-stellar cores (e.g., Bergin & Tafalla|[2007[ 
jCaselli & CeccareIIi||2012| ). [Guilloteau et al.j p006] l haw since 
detected DCO"^ in DM Tau, with a lower ratio of /?d ~ 4 x 10“^. 
A more recent, higher angular resolution study ofjQi et al.|(|2008]l 
has shown that the DCO^/HCO^ ratio in the TW Hya disk in¬ 
creases radially from 0.01 to 0.1 up to a radius of sa 9 0 au, 
where it drops off considerably. Later, Oberg et al. (2012 1 have 
observed isotopologues of HCN and HCO^ in TW Hya with 
both then SMA and ALMA, finding that the radial distribution 
of DCO^ and DCN is markedly different. While DCN seems to 
be centrally peaked, DCO^ shows an increasing column density 
with radius. This supports the theoretical predictions that these 
two deuterated ions are synthesized via distinct low- and high- 
temperature fractionation pathways. Recently, [Mathews et al. 


Another way to better characterize the thermal structure of 
PPDs is to observe and analyze deuterated species. Unfortu¬ 
nately, the key species for deuterium chemistry that can be de¬ 
tected at submillimeter wavelengths, namely ortho-H 2 D^ and 
para-D 2 H^, have been only observed in cold prestellar cores 


( 2013[ l have directly imaged the location of the CO snowline in a 
warmer disk around a Herbig A3 star HD 163296 with ALMA, 
using the optically thin DCO^ (5-4) line as a direct tracer of CO 
(see also [Qi et al.|2013| ). 

Additional information from the analysis of the HCO^ and 
DCO^ data includes the possibility to better c onstrain the io niza- 
tion degree than with HCO^ data alone (e.g., Caselli|2002 1 . The 
ionization degree is a key quantity that enables angular momen¬ 
tum transport in disks, thereby regulating their overall evolution 
and the ability to form planets via turbulence. Magnetorotational 
instability (MRI; Balbus & Hawley ( 1991[l) is currently the most 


widely accepted source of turbulence in disks (see, e.g. Flock 
[et al.|20T2|), although magnetocentrifugal disk winds and other 


in disks other more readily observable tracers, such as DCO 
DCN, DNC, N 2 D^, along with their major isotopologues, must 
be used. 

We must use a combination of these molecules as there are 
two main deuteration pathways possible, each with a different 
range of temperatures where they are most efficient. For exam¬ 
ple, DCO^ and N 2 D^ fractionation occurs mainly via H^ iso¬ 
topologues at temperatures < 20- 30 K, whereas fractionation 
of DCN and DNC involves deuterated light hydrocarbon ions 
such as CH 2 D^ and C 2 HD^ with a pathway which remains ac¬ 
tive up to temperatures of 70 - 80 K (e.g., [Millar et'ar]|1989[ 


non-lin ear effects such as am bipolar diffusion can also be im¬ 
portant. Oberg et al. (201 lb I have used the CO, HCO^, DCO^, 
and N 2 H^ lines observed in the disk of DM Tau with SMA and 
estimated ionization degree through its molecular layer. In the 
region probed by HCO^ (T > 20 K) the ionization degree was 
found to be 4 X 10 whereas in colder, deeper layers where 
N 2 H^ and DCO^ abundances are peaked, the ionization degree 
is lower, ~ 3 x 10“". 

The goal of this paper is to build on the previous work of 


Oberg et al. (201 Ibj) to better understand the thermal and ion¬ 


ization structure of DM Tau. We will present higher resolution 
observations of HCO+ J = (3-2), J = (1-0) and DCO+ J = (3-2) 
which, combined with a chemical model of DM Tau, allow us 
to determine the radial dependance of HCO^ and DCO^ column 
densities. Furthermore, the combination of two rotational lines 
for HCO^ with different optical depths provide us with a tool 
with which to discern the HCO^ excitation structure of the disk. 
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The paper is structured as follows: Sect. will describe the 
treatment of the observational data and the method in deriving 
a model. Sect. [^details the computational model used and a de¬ 
scription of the resulting ‘best-fit’ model to our observations. In 
Sect. 1^ we explore the deuterium fractionation and ionization 
fraction in the disk with a suite of chemical models to aid analy¬ 
sis, before finally summarising our results in Sect. 


2. Observational Results 

This section describes the process of creating a model of DM Tau 
by fitting HCO^ and DCO^ line emission. 


2.1. PdBl Data 


Observations were carried out with the IRAM Plateau de Bure 
interferometer. Table [TJpresents basic stellar properties and disk 
parameters of DM Tau from previous studies. We observed two 
transitions of HCO^, the J = (1-0) line at 89.18852 GHz and the 
J = (3-2) line at 267.55762 GHz and the J = (3-2) transition of 
DCO+ at 221.611258 GHz. The HCO+ J = (1-0) data included 


those described in Pietu et al. (2007 1 and were completed by 
longer baselines data (baselines up to 760 m) obtained in late Feb 
and early March 2008. The DCO^ data were obtained between 
Aug and Dec 2007, with baselines ranging from 15 to 175 m, 
yielding an angular resolution around 1.3". The HCO^ J = (3-2) 
data were obtained in Dec 2008, with a similar baseline coverage 
than the DCO^ data. 


2.2. Data Reduction 

We used the IRAM package GILDA^ for data reduction and 
imaging. All data were smoothed to similar spectral resolutions, 
0.17 to 0.20 km s“' for best comparison. Self-calibration was 
applied to all three lines and the dust thermal continuum was 
subtracted from the line spectra ( Dutrey et al.|2007' I. 

Channel maps of the emission lines are shown in Fig.[T] with 
contours in levels of 2.5cr for the HCO^ data and 2cr for DCO^. 
The rms noise values were calculated in a line free channel and 
are 3.7, 100 and 30 mJy/beam respectively. The need to plot 
tighter 2 cr contours for DCO^ J = (3-2) highlights the lower 
intensity of this line relative to the HCO^ emission, which is to 
be expected for a deuterated isotopologue. 


2.3. DiskFit Fitting 

Following the prescription of |Pietu et al.| ( |2007] ), we assume that 
the physical properties which affect line emission from a disk 
vary as a radial power law: 


a(0 = ao(^) . (1) 

where ao is the parameter value at the reference radius Rq- We 
adopt the standard that positive exponents, ea, imply a decrease 
of the physical quantity with radius. Through this prescription 
line emission in a disk can be described by the following param¬ 
eters: 

- Ao (") and Yq ("), position of the central star; 

- Vlsr (kms^'), systemic velocity; 

^ http://www.iram.fr/IRAMFR/GILDAS 


- PA (°) and i (°), position angle and inclination of the disk; 

- y„ (kms^'), (au), v, rotation velocity at reference radius 

R„, typically 100 au, with a power law exponent v. Perfect 
Keplerian rotation would yield v - 0.5; 

- Tm (K), Rt (au), qm, gas temperature at the reference radius 
Rj and respective power law exponent; 

- dy (kms '), e„, turbulent component of the line widths and 
respective power law exponent; 

- 2„, (cm“^), R^ (au), p^, gas surface density, reference radius 
and power law exponent; 

- Rin and Rout (au), inner and outer radius of observed emis¬ 
sion; 

- hn, (au), R/i (au), e/,, scale height of the gas, reference radius 
and its exponent. 


Note that the inclination and position angle are chosen in the 
range 0 < PA < 360 and -90° < i < 90° such that ymo is 
always positive and that position angle refers to the rotation axis. 

It has been shown previously that power laws are reasonably 
accurate proxies for correct descriptions of the large-scale disk 
properties, such as kinematics (assuming that the self-gravity 
is negligible), kinetic temperature distribution, surface density 
(assuming the a-viscosity prescription and a constant accretion 


rate), and thus, also the disk scale height (e.g., Chiang & Gol- 
|dreich|1997[|Pietu et al.|2007] i. 

The corresponding disk parameters were fitted to each of the 
obser ved lines with the DISKFIT software ( |Guilloteau & Dutrey| 


1998 1 , using a combination ofx^ minimization and MCMC fit¬ 


ting of the observed visibilities in the uy-plane. Line emission 
provided information on the gas structure while continuum emis¬ 
sion was used to fit the dust structure. 


2.4. Observationai Resuits 


The derived best-fit parameters for the three emission lines and 
the continuum are shown in Table Geometrical properties of 
the disk (systemic velocity, inclination, position angle and Kep¬ 
lerian velocity at 100 AU) were found to be in good agreement 
with previous studies on DM Tau using the DISKFIT approach 


(Guilloteau & Dutrey 1998[ Pietu et al. 2007 [ Dutrey et al.|2007| l. 


Due to the uncertainty in optical depth of the HCO^ rota- 
tional lines, we analyzed the HCO^ images in three different 
ways. In case A, the J = (1-0) and (3-2) images were fitted in¬ 
dependently. The derived parameters are presented in Col. 1-2 
of Table The HCO^ J = (3-2) is largely optically thick, and 
thus provides a good estimate of the temperature, however, the 
surface density cannot be well constrained from this data. A best 
fit outer radius of around 500 AU was found for this line. On the 
other hand, HCO^ J = (1-0) extends further (to at least 750 AU), 
and is mostly optically thin. Thus, the derived temperature heav¬ 
ily relies on the power law extrapolation at low radii. We also 
note that the J = (3-2) line requires an inner radius around 50 AU. 
A lack of emission from the inner 50 AU may bias the tempera¬ 
ture derived from the J = (1-0) line towards low values. 

In case B, we assume the J = (1-0) and (3-2) lines to have 
the same excitation temperature. We fit both transitions simul¬ 
taneously, and set the outer radii to 750 AU. Results are given 
in Col. 3. The derived temperature law is now much steeper, 
q - 1.00 ± 0.04, and the surface density law much flatter. The 
derived surface densities are also lower roughly a factor of 2 to 
5 between 200 and 400 AU compared to the separate fits. 

Finally, in case C, we relaxed the power law assumption, and 
fitted the temperature and surface densities at 5 different radii, 
extrapolating by power laws in between. The solution is within 
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4 ^ 33 '" 49!‘0 48"5 


Right Ascension 


Fig. 1: Channel maps of; HCO^ J = (1-0), top, with cr - 3.7mJy/beam (0.32 K) and contour spacing of 2.5cr; HCO^ J = (3-2), 
middle, with cr = lOOmJy/beam (0.64 K) and contour spacing of 2.5cr; and DCO^ J = 3-2, bottom, with cr = 30mJy/beam 
(0.42 K) and contour spacing of 2cr emission from DM Tau. The beam size, shown by the hlled ellipse, for each line is 1.4 x 1.26", 
1.85 X 1.49" and 1.4 x 1.33" respectively. The velocity of each channel shown in the top left corner in kms“^ while the cross in 
the centre shows the position of the host star and the position angle of the major and minor axes of the disk. Dashed contours show 
negative values. 


the errors consistent with case B, and thus we only report this 
latter case hereafter. 

Neither approach is perfect: in case A, the extrapolation re¬ 
quired to derive the temperature prohle from the J = (1-0) line is 
hazardous. On the other hand, in case B we neglect the vertical 
temperature gradients which are expected in disks, despite the 
higher opacity of the J = (3-2) line naturally leading to higher 
excitation temperature compared to the J = (1-0) line. Alterna¬ 
tively, if the density is insufficient to thermalize the J = (3-2) 
transition, we may expect its excitation temperature to be lower 
than that of the J = (1-0). We shall use the results of case B here¬ 
after. The very low temperatures derived in the outer part suggest 
some sub-thermal excitation, at least for the J = (3-2) transition 
beyond 400 AU or so. Note that apparent “ring-like” distribution 
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of HCO^ J = (3-2) is not due to simple excitation effect, as sug¬ 
gested by |Cleeves et al.| ( |2014| ). A central hole of 50 AU radius, 
almost fully devoid of HCOu is required to reproduce the J = 
(3-2) emission. 

The DCO^ line is weaker, and we can only derive the surface 
density by fixing all other parameters guided by those found for 
HCO^ as it is a reasonable assumption to assume they are cospa- 
tial in the disk. Note, however, that the derived values are quite 
insensitive to the assumed temperature law. 

To provide a better comparison with the modeled results, we 
will extrapolate the power laws describing column densities and 
temperatures. However it must be noted that, due to the inner 
hole in HCO^ and smaller outer radius in DCO^ the power laws 
only provide a good fit to the data in the region 50 ^ r ^ 430 AU. 
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Table 2: Best-fit parameters for DM Tau with descriptions of the parameters in the text. 


Parameters 

HCO+ J = (1-0) 

HCO+ J = (3-2) 

HCO^ Simultaneously 

DCO+ J = (3-2) 

Continuum 

Vlsr 

(kms *) 

6.05 ± 0.01 

6.01 ± 0.02 

[6.01] 

6.00 ±0.21 

[6.01] 

i 

(°) 

34.0 ± 2.7 

33.8 ±0.5 

[34] 

34.5 ±2.1 

[34] 

PA 

n 

64.31 ±0.57 

65.7 ± 3 

[65] 

65.9 ± 1.3 

[65] 

Vioo 

(kms^*) 

2.06 ±0.10 

2.16 ±0.05 

[2.1] 

[2.1] 

[2.1] 

^int 

(an) 

- 

53 ±7 

49/^ 

70 ±20 

- 

^out 

(an) 

[750] 

510±5 

[800] 

427 ± 10 

173.5 ± 0.3 

dV 

(kms^*) 

0.17 ±0.01 

0.14 ±0.02 

- 

0.22 ± 0.44 

[0.15] 

him 

(an) 

[16.5] 

[16.5] 

[16.5] 

[16.5] 

[16.5] 

Tim 

(K) 

11.3±0.2 

19.0 ±0.2 

33.6/it 

[17] 

19.07 ±0.04 

1 


0.34 ± 0.04 

0.46 ± 0.04 

1.00 ±0.04 

[0.43] 

0.44 ± 0.03 

O 

O 

(cm“^) 

(1.9 ± 0.8) x 10*^ 

[ 2.00 x 10*^1 

( 9 . 8 !“) x 10*2 

(1.2 ± 0.7) x 10*2 

(2.2 ± 0.66) x 1022 

P 


2.6 ± 0.8 

[2.5] 

0.82 ± 0.06 

0.44 ±0.11 

0.61 ±0.05 


Notes. Rotational lines result in parameters for the gas, the continuum probes the dust. Values in square brackets were fixed during fitting. 



40 100 600 

Radius (AU) 

Fig. 2: Disk physical structure scaled vertically as zjR. Colour¬ 
ing shows the coupled gas and dust temperature and solid lines 
show the particle number densities in particles cm“^. The dashed 
white lines shows the 21 K isotherm, below which half of CO 
has frozen out. 


In order to better comprehend these findings, we perform de¬ 
tailed theoretical modeling in the next Section. 


3. Computational Model 

3. 1. Disk Physical Structure 


This section describes the methodology of creating a computa¬ 
tional chemical disk model of DM Tau, including a description 
of the physical and chemical parameters used. 

The DM Tau system, at a distance of 140 pc, consists of a sin¬ 
gle isolated pre-main-sequence MO.5-1.5 dwarf (Tefr = 3720 K), 
with a mass of 0.5 - 0.65 Mq, a r adius of 1.2 R^, and an accre- 
tionrateof~ 2-3x 10“® Mq yr“' (Mazzitelli 1989 Simonetal. 


2000[ McJunkin et al.p014i. It is enshrouded by an extended 


(~ 800 AU), cold (T ^ 10 K) Keplerian disk (|Pietu et al. 


According to the Spitzer IRS observations (Calvet et al. 


2007 


2005 


the inner DM Tau disk is cleared of small dust (;^ 3 - 4 AU) 
and is in a pre-transitional phase. As our interferometric obser¬ 


vations have the highest sensitivity in disk regions > 30 AU from 
the central star, we only consider the chemical evolution outside 
of this radius in our analysis. 

The DM Tau physical disk model is based on a 1-i-lD steady- 
state cr-model similar to that of |D’Alessio et al?] ( |1999| ), where 
equal gas and dust temperatures are assumed. This model was 
extensively used in our previous studies of DM Tau-like disk 
chemistry (e.g., Henning et al. 2010 Semenov & Wiebe |201 1 


[Albertsson et al.|2014b| l. The disk model has an outer radius of 

800 AU, an accretion rate of 2 x lO^^Mgyr^^ a viscosity pa¬ 
rameter a - 0.01, and a total gas mass of 0.066 Mq ([Dutrey 
|et al. |2007[ [Henning et al.|[2()10[ Semenov & Wiebe||2011| l. 
The dissociating UV radiation of DM Tau is represented by the 


scaled-up interstellar UV radiation field of Draine (19781. The 


unattenuated stellar UV intensity at the radius of 100 AU is 
;^^t(100) = 410 (e.g., [Bergin et al.]|2003| l. The X- ray luminos 


ity of DM Tau is taken to be 2 x 10" ergs ‘ (see Semenov & 


|Wiebe|201 Ij l. The calculated disk thermal and density structure 
is shown in Fig.[^ 


3.2. Disk Chemical Model 

The adopted chemical model is based on the advanced 
ALCHEMIC code (see Semenov et al.|20f0| ) and utilizes the high- 
temperature, gas-grain deuterium chemistry network of Alberts 


son et al. ( 2013[ ), with the addition of nucle ar spin-state pro- 
and H; 


cesses for Ht, Hf 


Albertsson et al. 


12 , ti 2 , anu rij isotopologues from 
(|2014al. The chemical network without deuterated species 


based on the osu.2007 ratefil^ with the recent updates to the re¬ 
action rates from Kinetic Database for Astrochemistry (KIDA) 
( jWakelam et al. |2012| ). For all H-bearing reactions in this net¬ 
work, we derived the corresponding D-bearing reactions follow¬ 
ing the algorithm of [Rodgers & MilE| ( |1996| . The cloning was 
not allowed for any species with the -OH endgroup. 

Primal isotope exchange reactions for H:t as well as CHJ 
and CzH^ ~ ~ ' 


from 


_ Roberts & Millar (20001; Gerlich et al. ( 2002|; 

[Roberts et al.| ( |2(J04| ); |Kouett et al.| ( |2005| l vvere included. In cases 
where the position of the deuterium atom in a reactant or in a 
product was ambiguous, a statistical branching approach was 


See: 

research.html 


http: //WWW. physics. Ohio- state. edu/~eric/ 
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INITIAL ABUNDANCES FOR MODELING PRE-DISK EVOLUTIONARY PHASE 


ortho-H 2 

0.375 

He 

9.75(-2) 

0 

1.80(-4) 

Na 

2.25(-9) P 

2.16(-10) 

para-H2 

0.125 

C 

7.86(-5) 

S 

9.14(-8) 

Mg 

1.09(-8) Cl 

1.00(-9) 

HD 

1.55(-5) 

N 

2.47(-5) 

Si 

9.74(-9) 

Fe 

2.74(-9) 



25 INITIALLY MOST ABUNDANT MOLECULES FOR DISK CHEMICAL MODELING INCLUDING HCO^ AND DCO^ 


para-H2 

3.77(-01) 

CO* 

4.05(-05) 

nh; 

5.64(-06) 

H* 

6.03(-07) 

NO 

2.22(-07) HCO+ 

6.13(-09) 

ortho-H2 

1.23(-01) 

CO 

3.26(-05) 

0 

5.59(-06) 

CjH; 

4.48(-07) 

N 

1.36(-07) DCO+ 

1.25(-11) 

He 

9.75(-02) 

02 

1.79(-05) 

0 ; 

4.12(-06) 

OH 

3.43(-07) 

HDO* 

1.35(-07) 


H 

5.25(-04) 

HD 

1.52(-05) 

ch; 

3.64(-06) 

H 20 

2.79(-07) 

n 

p 

1.32(-07) 


H 2 O* 

5.53(-05) 

N 2 

7.39(-06) 

N* 

1.76(-06) 

HNO* 

2.40(-07) 

n 

P 

1.19(-07) 



Table 3: Atomic and molecular abundances used in the disk chemical model. Top; Initial abundances for modeling the pre-disk 
evolutionary phase Bottom: 25 initially most abundant molecules for disk chemical modeling. Note that a{b) should be read as 
a X 10* and that * denote frozen species. 


used. This deuterium network was further extended by adding 
ortho- and para-forms of H 2 , and isotopologues and the 
related nuclear spin-state exchange processes from several ex 


perimental and theoretical studies (Gerlich 

1990 

Gerlich et al. 

2002 

[Flower et al.||2004[ [Walmsley et al. 

2004 

Flower et al. 

2006 

Pagani et al.[2009[ Hugo et al.[2009 

Tonvault et al.[201 1 

Sipilaetal.[2013[l. 


To calculate UV ionization and dissociation rates, the mean 
FUV intensity at a given disk location is obtained by summing up 
the stellar= 410x(r/100)“^, where r is in au, and interstel¬ 
lar UV fluxes scaled down by the visual extinction in the radial 
and vertical directions, respective ly. Several tens of newer ph o- 
toreaction rates are adopted from van Dishoeck et al. (2006 1 ^ 


The self-shielding of H 2 from photodissociation is calculated by 
Eq. (37) from Draine & Bertoldi| ( |1996| l. The shielding of CO by 


dust grains, H 2 , and the CO self-shielding is calculated using a 
precomputed table of |Lee et al.| ( fT99^ Table 11). 

The stellar X-ray radiation is modeled using observational 
results of Glassgold et al. (2005|l and the approximate expres¬ 


sions (7-9) from the 2D Monte Carlo simulations of Glassgold 


|et al. (|1997a|b|l. Implementing Eqn. (8) from [Glassgold et al. 


1997b I, we use an exponent of n = 2.81, a cross section at 
keV of (T -22 - 0.85 x 10“^^ cm^ and total X-ray luminosity of 
Lxr = 3 X 10^® erg s“', yielding a typical X-ray photon energy 
of 3 keV. Attenuation of X-rays is calculated from Eqn. (4) in 
[Glassgold et al.|p997a) l. The X-ray emitting source is located at 
12 stellar radii above the midplane with rates exceeding that of 
the CRPs in the disk regions above the midplane, particularly, at 
radii ~ 100 - 200 AU (see also [Henning et al.|2010[ l 

We assume the standard cosmic ray (CR) ionization rate 
^CR = 1.3 X 10“'^ s“' and model its attenuation using Eq. (3) 
from [Semenov et al.| ( [2004| l. Note that we do not consider 
the scattering of low energy CR protons by the heliosphere of 
DM Tau, as done in [Cleeves et al.| ( [2013a| l. Ionization due to 
the decay of short-livi ng radionuclides is ta ken into account, 
^RN = 6.5 X 10 s ' (Finocchi et al.||l997ll, see also Cleeves 
[etakldlOT^ . 

The grain ensemble used to calculate disk physical struc¬ 
ture was represented by uniform amorphous silicate particles 
of olivine stoichiometry with density of 3 gcm“^ and ra dius 
of 0.1 pm. Each grain provides a; 1.88 x 10* surface sites (Bi- 
[ham et al.|2001[ l for surface recombination that proceeds solely 
through the classical Langmuir-Hinshelwood mechanism (e.g. 


http: //WWW.strw.leidenuniv.nl/~ewine/photo 
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[Hasegawa et al.|1992| l. The gas-grain interactions include stick¬ 
ing of neutral species and electrons to dust grains with 100% 
probability and desorption of ices by thermal, CRP-, and UV- 
driven processes. We do not allow H 2 to stick to grains as it re¬ 
quires temperatures of $ 4 K. The UV photodesorption yield of 
3 X 10“^ was adopted (e.g., Oberg et al.|2009a b} Eayolle et al. 


[201 1[ [2013|l. Photodissociat ion processes of solid species are 
taken ffom [Garrod & Herbst| ( [20()^ ; |Semenov & Wiebe|P0lT] l. 

In addition, dissociative recombination and radiative neu¬ 
tralization of molecular ions on charged grains and grain re¬ 
charging are taken into account. Upon a surface recombination, 
we assume there is a 1% probability for the products to leave 
the grain due to the partial conversion of the reaction exother- 


micity into breaking the surface-adsorbate bond (Garrod et al. 
2007[ Vasyunin & Herb^[2013[l. Eollowing experimental stud¬ 


ies on the formation of molecular hydrogen on amorphous dust 
grains by Katz et al. ( [1999[ ), the standard rate equation approach 
to the surface chemistry was utilized. Overall, the disk chemi¬ 
cal network consists of 1268 species made of 13 elements and 
38812 reactions. 

The age of the DM Tau system is poorly constrained, ~ 3- 
7 Myr ( [Simon et al. [2000] ). In the chemical modeling the age 
of 5 Myr was considered. To set initial abundances, we calcu¬ 
lated chemical evolution in a TMCl-like molecular cloud {tin = 
2x10"' cm“^, T = lOK, Av = 10 mag) over 1 Myr. Eor that, 
the “low metals” elemental abundances o f|Graedel et al.| ( |l982[ l; 
[Lee et ah ( [1998| l; [Agilndez & Wakelam ( [2013[ l were used, with 
the equilibrium 3:1 ortho/para H 2 ratio (hydrogen being fully in 
molecular form) and deuterium locked in HD molecule (see Ta¬ 
ble]^. The resulting abundances of modeling the pre-disk evolu¬ 
tionary phase were used as initial abundances for disk chemical 
modeling as shown in Table 


3.3. Modelled Results 

Eigure.j^ presents the modeled column densities for HCO^ and 
DCO^ in DM Tau (dashed lines) and associated 3 cr errors with 
the observationally derived values (solid lines) overlain. Both 
column densities agree well within their errors. The grey boxes 
in Eig[^show where the power-law column densities are extrap¬ 
olated beyond the radii where they were observed as to provide 
a better comparison with the chemical model which does not re¬ 
produce the inner hole. 

Eigures and c show the relative abundance with respect 
to H 2 of DCO^ and HCO^ throughout the disk. These clearly 
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Radius [AU] 

Fig. 3: Deuterium fractionation in the disk, (a) Comparisons of HCO^, blue, and DCO^, red, from observations, dashed with 3cr 
errors, and the best fit model, solid with associated error of a factor of 3. They grey vertical bars show where the observationally 
derived column densities are extrapolated beyond the inner and outer edges found in parameterisation. (b) and (c) show the relative 
abundances of DCO^ and HCO^ in the best fit disk chemical model respectively, (d) shows the local DCO^/ HCO^ ratio in the disk 
model. 


demonstrate the stratification of the disk with a distinct molecu¬ 
lar layer lying ~ 0.5 pressure scale heights above the midplane. 
Both molecular layers are relatively co-spatial, in general trac¬ 
ing regions of high gas phase CO abundance. They are bounded 
by the CO snowline towards the midplane, and the photodisso¬ 
ciation region of CO towards the disk atmosphere. DCO^ occu¬ 
pies a slightly tighter vertical range than HCO^ due to it also 
requiring efficient deuteration. The molecular layer is truncated 
upwards due to the higher gas temperatures, reducing the effi¬ 
ciency of deuterium fractionation, while the lower bound is due 
to reduced HD abundances with which to readily transform deu¬ 
terium into DCO^. 

/?d(HCO^) is shown in Fig. w with local values reaching 
as high as ~ 1. Within this region we see pronounced increases 
of HD and H 2 D^ abundances facilitating a fast transfer of deu¬ 
terium from HD to DCO^. Despite these locally high values, 
the vertically integrated column density is more sensitive to the 
denser regions closer to the midplane, thereby exhibiting a lower 
value of/?D(HCO+)~ 0.1. 

Figure shows the modelled values of Ro a various time 
steps in the model (blue lines). The observationally derived 
value, shown by the black dashed line, agrees qualitatively well 
with the f = 5 X 10® yrs value. 


4. Discussion 

In this section we use both the observationally-derived and best- 
fit modeled column densities, complemented with a further suite 
of chemical models, to explore both deuterium fractionation and 
the ionization fraction in the disk. 


DCO+ Molecular Layer Temperature (K) 

25 20 15 10 



Fig. 4; Comparisons of /?d(HCO^) values. The dashed black 
line shows the observationally derived value in DM Tau and 
the blue solid lines show the chemical model values at differ¬ 
ent time steps. Typical errors are a factor of 3. They grey verti¬ 
cal bars show where the observationally derived column densi¬ 
ties of HCO^ and DCO^ are extrapolated beyond the inner and 
outer edges found in parameterisation. The orange dash-dotted 
line shows the value observed in TW Hya dQi et al.||20081l, the 
red dotted line, /?d(HCO+) in HD 163296 ( [Mathews et al.|2013| l, 
and the yellow dashed region, values from a survey of prestellar 
cores ( [Butner et al.|1995} . The top y-axis shows the temperature 
of the DCO^ layer derived from our parametric fitting. 


4.1. Deuterium Fractionation of HCO^ 

Deuterium fractionation is typically used as a probe of the ther¬ 
mal history of an environment. In this section we explore phys¬ 
ical parameters other than temperature which can alter the ob¬ 
served Ro = Al(DCO^)/ A^(HCO^) ratios in disks through a suite 
of chemical models derived from our DM Tau model. 

Figure ^shows the radial profile of Ro found for DM Tau ob¬ 
servationally (dashed black line) and at various time steps in the 


model, t - {10^, 10"^, 10®, 10®, 5 x 10®) yrs (solid blue lines). 
Also marked are values for TW Hya (orange dash-dotted line; 
IQi et al.| p008| l), HD 163296 (red dotted line; [Mathews et dT] 
( 2013|l) and an av erage of prestellar cores (yellow hatched re¬ 
gion; Butner et al.[ ( [1995] l; [Caselli et al.[ ( [2002[ l). Enhanced values 
relative to the cosmic abundance of [D]/[H] ~ 10 ® are indica¬ 
tive of continued gaseous processing in a cold, T ^ 20 - 30 K, 
environment; a none too surprising conclusion given the kinetic 


Article number, page 7 of 


13 















































A&A proofs: manuscript no. hco_dco_obs 



R, AU 

(a) Increased grain size. 



R, AU 

(b) Increased X-Ray luminosity. 


R, AU 

(c) No stellar UV. 



100 
R, AU 

(d) No CRPs. 



R, AU 

(e) Reduced X-Ray Luminosity 


R, AU 

(f) No interstellar UV. 


Fig. 6: Computed molecular column densities for DCO^ (dashed line) and HCO^ (solid line). Black lines show the columns obtained 
with the best-fit model (see Section 3.3 i, red lines show the columns when the specified physical parameter is changed in the best 
ht model while the others are held constant. For reference, the canonical model has a single grain population of a = 0.1 pm. 
Lx = 3 X 10^® erg s-',;i^*(100 AU) = 410 and a CRP ionization rate of ^ = 1.3 x 10^'^ s^'. 



40 100 600 

Radius [AU] 


Fig. 5: (a) The A(ortho-H 2 )/A(para-H 2 ) ratio from the best ht 
model. After ~ 10^ yrs the ratio reaches a steady state, (b) 
The degree of CO freeze out in our model, A(COice)/A^(COtotai) 
clearly demonstrating the CO ice line at r ^ 30 AU. 


temperature probed by the line emission, T ~ 10 - 20 K (see 
Table ig. Furthermore, the radial increase is to be expected due 
to the radial temperature gradient in a disk. That is, the pace of 
deuterium fractionation and synthesis of the isotopologues 
hastens in colder outer disk regions. In turn, higher abundances 
of the isotopologues imply more efficient formation of DCO^ 
in ion-molecule reactions with gaseous CO, increasing Rq. 


Our Ru values, ranging from « 0.1 - 0.2 between 50 AU 
and 430 AU, are almost two orders of magnitude higher than 
the disk average value of (4.0 + 0.9) x 10^^ found by 


Guilloteau 


et al. (20061, who used HCO J = (1-0) interferomic data and 
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DCO^ J = (3-2) single dish data. This discrepancy can be ex¬ 
plained by the assumptions made in the calculation of HCO^ 
and DCO^ column densities. Firstly, it was assumed that DCO^ 
was radially co-spatial with the HCO^ emission which extends 
out to ~ 800 au, an outer radii similar to that of CO. This in¬ 
creased emitting region would result in a lower, disk average 
value for DCO^, reducing R^. Secondly, we have shown that 
HCO^ likely exhibits a complex molecular distribution includ¬ 
ing an inner hole in emission, see Sect. |2.4| When fitting a single 
transition this complexity was found to drive the parameterisa- 
tion to favour steeper power laws describi ng column densities , 
consistent the HCO^ J = (1-0) data from Pietu et al. ( 2007 1 , 
yielding a much smaller Ru value in the inner disk than found 
with our data. 

In TW Hydrae, another well studied protoplanetary disk, 
the value of Ru has been measured both as a disk average 


(0.035 + 0.015, van Dishoeck et al. (2003l), and with spatially 
resolved interferomic observations yielding a value from 0.01 to 
0.1 between 30 and 90 au ( |Qi et al. 12008) 1. The increase found in 
the outer regions of the disk when moving to spatially resolved 
data is consistent with our findings in DM Tau. Similarly, while 
the TW Hya disk is smaller and less massive than DM Tau, it 
appears to hold comparable values for Ru at similar distances to 
the star, namely around ~ 0.1 at 100 AU. More recently, |Math-| 
lews et ah] ( |2013| ) used ALMA science verification data of HD 
163296, an Al spectral type Herbig Ae star and comprehensive 
modeling to ascertain local values of Ru ~ 0.3 and a disk aver¬ 
age of 0.02. Typically disks around Herbig Ae stars are warmer 
and more massive than those around classical T-Tauri stars, thus 
a reduced disk average value of Ru is to be expected. 

Figurej^also clearly shows continued enhancement of Ru up 
to ~ 1 Myr, indicative of continued processing of the gaseous 
CO during a disk lifetime. Thus, it is of no surprise that the 
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observed value of in DM Tau, thought to be between 3 - 
7 Myr, is higher than that found in low-mass prestellar cores 
(0.045 + 0.015, |Butner~erar] ( fT995l )). 

Molecular layers in a protoplanetary disk are far 
beyond closed systems. Measurements of Rq through 
A^(DCO^)/A^(HCO^) are sensitive to more than just the 
deuterium fractionation efficiency, but are compounded by 
physical parameters which can alter the abundances of HCO^ 
and DCO^. In the remainder of this section we explore how 
else the physical environment can change the observed value of 
Rd by running a suite of chemical models and varying a single 
parameter. This will aid in analysis when comparing values of 
Ru from different astrophysical environments. To motivate the 
choice of parameters studied, we discuss briefly the process of 
deuterium fractionation. 

Deuterium fractionation occurs due to the energy difference 
between deuterium and hydrogen atoms, thereby resulting in the 
deuteration of Hj by HD to be exothermic: 


HD + H7 H2 -H H 2 D+ -H AE, 


( 2 ) 


where AE - 232 K ( [Roberts & Millar|2000| l. Hence, low energy 
environments are conducive to enhanced abundaces of deuter- 
ated isotopologues resulting from successive deuteration of H^. 
For H 2 , AE corresponds to Tkin ~ 30 K resulting in inefficient 
fractionation above this temperature. This is clearly evident in 
both observational and modelled results (see Fig.[^ which show 
an increase in Rd at larger radii where the disk is cooler. Further¬ 
more, the larger internal energy of ortho-H 2 allows for greater 
leverage in the back reaction of Eqn. It has been shown that 
ortho / para ratios >0.1 can ‘poison’ the overall fractionation 
efficiency, even when the kinetic temperatures are low. Proton 
exchange due to collisions will reduce the initial H 2 ortho / para 
fraction of ~ 0.75. A reduced ortho-para fraction will increase 
the efficiency of deuterium fractionation, see Fig|5^. (e.g.. Pa 


gani et al.|1992||2009[ Crabtree et al.|2011[ Pagani et al.|2()13[ 


Albertsson et al.| 2014al. However, the relatively high densities 

of the molecular layer in a disk ensure that this is a relative quick 
process taking ~ 10^ yrs. 


4.1.1. Regulating CO Abundance 

As a parental molecule of HCO^ and DCO^, the CO abundance 
is intimately linked to the abundances of HCO^ and DCO^. 
Gaseous CO must be sufficiently abundant to efficiently convert 
the Hj and H 2 D^ into HCO^ and DCO^ respectively. However, 
CO readily freezes out at T ^ 21 K, vastly reducing the avail¬ 
able reaction partners, see Fig|^. Therefore, DCO^ is most ef¬ 
ficiently produced where temperatures are high enough to main¬ 
tain a relatively low level of CO depletion, yet cool enough to 
allow efficient fractionation. This is visible in Fig. where 
A(DCO^) peaks around 50 AU. 

For the same mass of dust, larger grains have a reduced sur¬ 
face area onto which CO can freeze out, reducing the deple¬ 
tion of CO and hence expanding its molecular layer towards 
the disk midplane. Fig. 6a shows the change in A(HCO^) and 
A(DCO^) when the grain size in our best fit model is increased 
to 1 jum. Outside the CO snowline A(HCO^) is enhanced due 
to the greater availability of CO to react with Hj^ resulting in a 
shallower gradient. On the other hand, A(DCO^) is uniformly 
decreased. The decrease of DCO^ abundances is associated with 
less efficient fractionation of H^ in its molecular layer. Rapid 
ion-molecule reactions of H7 with volatile molecules such as 


CO, which are less depleted from the gas phase, compete with 
fractionation processes and hence lower production of deuter- 
ated Hj isotopologues. This results in a reduced Rd as a disk 
average and a weaker radial gradient. 


4.1.2. Ionization 

The other parental molecules of HCO^ and DCO^ are H^ and 
H 2 D^ respectively, both of which require the ionization of H 2 . 
Figs |6bP| show how the column densities are affected when ion¬ 
ization sources in the best fit model are altered. 


Clearly shown in Figs. 6c androdl stellar UV and CRPs play 


little role in the abundance of these two species. In our disk 
model, UV scattering is neglected therefore the stellar UV radi¬ 
ation becomes quickly absorbed in radial direction by the dust, 
thus making little impact on the abundance of HCO^ and DCO^. 
Additionally, high energy CRPs can penetrate to the molecular 
layers, however they have such a low flux that they are of little 
consequence in the life-cycle of HCO^ and DCO^. 

Conversely, the abundances of HCO^ and DCO^ are sen¬ 
sitive to the stellar X-rays, the dominant ionization source in 
the molecular layer. Increased values of Lx lead to an enhanced 
abundance of HCO^ across the entire disk, as shown by Fig. 6b 
DCO^ production is suppressed in the inner disk, r < 100 AU, 
while outer regions display an enhancement. The DCO^ produc¬ 
tion is suppressed in the inner disk despite an increase in overall 
Hj isotopologue abundances within the inner molecular layer 
due to an increased ortho/para-H 2 fraction arising from the in¬ 
creased X-ray luminosity injecting sufficient energy into the disk 
for re-equilibriation through ion-molecule and nuclear spin-state 
processes. The increased ortho-fraction of the H^ isotopologues 
slows the overall pace of the deuterium fractionation. In addi¬ 
tion, an increase of DCO^ abundances in the outer disk is smaller 
than for HCO^ because the DCO^ layer is located deeper into 
the disk, where temperatures favor deuterium enrichment andare 
better shielded from impinging stellar X-ray photons than the 
more extended HCO^ molecular layer. 

Reducing the stellar X-ray luminosity reduced both the 
HCO^ and DCO^ abundances. This is a more pronounced ef¬ 
fect in the inner disk due to the oblique angles of incident X- 
rays, as found in Henning et al.| ( |20T0 i. The higher sensitivity 
of HCO^ to changing X-ray luminosities, again due to the dif¬ 
ference in vertical extents of the molecular layers, is reflected in 
the gradient of Rd; a lower Lx leads to a less radially dependent 
Rd. Additionally, HCO^ production is suppressed to such a low 
level that Rd can reach ~ 1 in the outer disk. 

A more puzzling result is the strong influence of interstellar 
UV radiation on Rd. In the disk model without an IS UV field, 
Rd becomes lower by up to an order of magnitude compared 
to the best-fit model. This is mainly due to a uniform decrease 
of the DCO^ column density throughout the disk by a factor of 
~ 3, and an increase of the HCO^ column density at r ^ 50 - 
60 AU by a similar factor of about 3, producing a near constant 
Rd across the radius of the disk. 

The IS UV photons play two main roles for disk chemical 
evolution. Firstly, they partly contribute to the ionization and 
dissociation of disk matter in the molecular layer, particularly 
beyond r ^ 100 AU. Secondly, and more importantly for Rd, 
they bring heavy ices back to the gas phase by photodesorption 
and thus partly regulate surface processes. We found that in the 
DM Tau model without IS UV radiation, CO gets more easily 
converted into CO 2 ice in the molecular layer at T 30 - 40 K 
through the slightly endothermic reaction of CO + OH + 80 K 
^ CO 2 + H. CO 2 is unable to be photodesorbed and therefore 
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partly dissociated in the gas. In addition, abundances of water 
ices increase, whereas abundances of atomic and molecular oxy¬ 
gen decrease. This leads to a drop in gas-phase CO abundances 
locally by a factor of 2-4 at all disk radii. 

Furthermore, in the absence of photodesorption due to the 
absence of UV ionization, ions of alkali metals such as Na^ and 
Mg^, along with atomic ions such as and C^, become less 
abundant and do not contribute considerably to the fractional 
ionization of the entire molecular layer. As a result, polyatomic 
ions like HCO^ and H 30 ^ dominate the ionization structure, re¬ 
sulting in a decrease by a factor of a few in ionziation fraction 
due to their more efficient recombination with electrons. These 
two factors, lower ionization degree and lower CO abundances, 
lead to both lower DCO^ abundances in the molecular layer and 
also reduced DCO^ column densities. 

Contrary to DCO^, HCO^ abundances and column densities 
show an increase at r ^ 50 AU in the DM Tau model with¬ 
out IS UV. This is related to the fact that the DCO^ molecular 
layer is narrower and located more deeply in the disk compared 
to the HCO^ molecular layer. The lack of photodesorption and 
photodissociation in the upper part of the HCO^ molecular layer 
increases abundances of CO and isotopologues, hence boost¬ 
ing production of HCO^. This compensates for the decrease of 
gaseous CO due to its surface conversion into the CO 2 ice in the 
lower part of the HCO^ molecular layer. 

Comparisons with previously published models highlight the 
importance of considering these additional processes which alter 
the abundance of HCO^ and DCO^. Despite the relatively well 
understood deuterium fractionation mechanism there are orders- 
of-magnitude disparity between models. [Aikawa et al.| P002| l 
modelled a smaller disk with gas extending on to 373 AU which 
exhibited a radial dependence of Rd varying between 0.003 to 
0.06, values more more in accord with the smaller TW Hya disk. 
Whereas a newer model of Willacy (2007 1 found considerably 


higher values ranging from Ru = 0.1 at 50 AU to reaching unity 
outside 100 AU, suggesting high Rq values are to be expected. 

Hence, while A(DCO^)/ A(HCO^) provides an easily acces¬ 
sible measure of deuterium fractionation, a link through several 
environments in the cycle of molecular gas, other parameters, 
particularly the X-ray luminosity of the central star, the inter¬ 
stellar UV field and the grain evolution, are folded into this mea¬ 
surement. In the case of protoplanetary disks, the X-ray lumi¬ 
nosity of the host star must be well constrained in order to fully 
characterise the deuterium fractionation present in the disk. 


4.2. Ionization Fraction, xie ) 

HCO^ is often touted as the most dominant ion in the warm 
molecular layer of a protoplanetary disk. As such, it is frequently 


used as a proxy of the ionization in this region dSemenov et al . 
2004| iDutrey et al.||2007| |Qi et al.|[2008l [Oberg et al.||2011b| l. 


However, the large radial and vertical gradients in physical pa¬ 
rameters characteristic of a protoplanetary disk introduce several 
complications in deriving knowledge of the ionization fraction 
from a single charged species. Common practice therefore is to 
make a steady state approximation, a methodology that has been 
applied to a range of astrophysical scales: from protoplanetary 
disks to molecular clouds and supernova remnants (|Guelin et al.| 


T977||Caselli|2002HCaselli et al |2008| [Vaupre et al.|20T^ 


Introduced by Guelin et al.| (|1977)l, this assumes a heavily 


reduced chemical network in a steady state of ionization. As dis¬ 
cussed previously, the abundance of HCO^ is largely governed 
by two main processes: creation through ion-neutral reactions 
between H^ and CO and destruction via electronic recombina- 
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Fig. 7: Constraints on the ionization fraction in DM Tau. Re¬ 
sults from the steady state approximation are shown with: obser¬ 
vational values, black dashed; modelled values from the entire 
disk, dark blue solid; modelled values from the HCO^ molecu¬ 
lar layer, light blue solid; lower limit frornjDuttey_e^alJ|2007 1 , 
light orange solid and lower limit from Oberg et al. (201 lb 1 , 
dark orange solid. All errors, dominated by the values from the 
chemical modelling, are a factor of 3. 
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D2H+ + CO ^ HCO+ + Dz 

kco 

1.61 X 10-® 

0 

- 

D2H+ -H CO ^ DCO+ -H HD 

kco 

1.61 X 10^® 

0 

- 

D+ -H CO ^ DCO+ + Dz 

kco 

1.61 X 10^® 

0 

- 

RECOMBINATION 



HJ -1- —» various 

^recO 

6.8 X 10^* 

-0.5 

- 

HzD^ + e —> various 

^recl 

6.0 X 

-0.5 

- 

DzH^ + e —» various 

^rec2 

6.0 X 10-® 

-0.5 

- 

Dj + —> various 

^recO 

2.7 X 10-* 

-0.5 

- 


Table 4: Reactions involved in the considered steady state system 
and the rates used in our modeling. Forward rates should be read 
as a(T/3QQ"f and backwards as a(T1300)^ exp{-ylT). Adapted 
fromICaselli et al.l(1200811 . 


tion. Similar pathwa ys hold for DC O 
HaD^ in place of H 3 . 
ics of such a network 


but with the deuterated 


Caselli (2002 1 showed the chemical kinet- 
m steady state can be reduced to: 


1 (kix(HD) \ 

x(e ) = - ——- kix(CO) , (3) 

«recl \ 3Ru ) 

where the rates and associated reactions are found in Tabled 
This requires knowledge of both V(HD), of which we only 
have one observation of in a protoplanetary disk, TW Hya 
(|Bergin et al.||201^, and the total gas column from which to 
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Fig. 8; (a) Ionization structure of the DM Tau disk model. The region bounded by blue lines is what we have defined as the HCO^ 
molecular layer, (b) and (c) Relative contributions of charged species to the ionization level over the whole disk and the HCO^ 
molecular layer. As a disk average, panel (b), it is clear the charge is dominated by and the gray region representing all 
other charged species. In the HCO^ molecular layer, panel (c), defined as the region that contributes 90% of the total HCO^ column 
density, HCO^ supplies a majority of the charge with large contributions from H^, and S^. 


convert from column density to relative abundances, a value 
that cannot be well constrained observationally without several 
assumptions (Thi et al. 2010| l. Application to DM Tau is thus 
limited to the case where we must assume values of x(CO) 
and x(HD). Note that this has since been expanded to include 
all multiply-deuterated isotopologues of and charged grains, 
however the lack of observational constraints on these would fur¬ 
ther compound the issues detailed above in the case of a proto¬ 
planetary disk ( Caselli et al.|2008' i. Note, however, that physical 
parameters derived from line emission will be of the molecular 
region and not applicable to the disk as a whole. 

Using our observationally derived Ru values and values 
x(HD) = 2.40 X 10“^ and x(CO) = 7.24 x 10“^ taken from the 
best fit model, the later consistent with previous observations of 
DM Tau (|Pietu et 'ar]|2007|l, we find an ionization fraction of 


x(e ) ~ 10 ', as shown by the dashed blac k line in Fig, [t] 
is consistent with the lower lirnits plac ed by|Dutrey et al.](|^ 


This 

^2007] l 


(light orange) and Oberg et al. (201 lb I (dark orange). Blue lines 
show x(e ) from our best fit model, the light blue considering 
molecular column densities integrated over the warm molecular 
layer probed by our HCO^ and DCO^ observation^ and those 
integrated over the whole disk in dark blue. Both of which quali¬ 
tatively agree with the steady state values. IQi et al.| ( |2008| l found 
a similar ionization fraction of x(e ) ~ 10 ' in TW Hya when 
using the same steady state approximation. 

However, disk ionization, is controlled by a myriad of atomic 
and molecular ions as shown by Fig[^ Panels (b) and (c) show 
the relative contribution of the top eight most abundant ions as a 
function of radius for the whole disk (b) and the HCO^ molec¬ 
ular layer (c). It is clear that ionization as a whole is dominated 
by the atomic ions and which contribute > 99 % of the 
charge. Even within the molecular layer, atomic ions are the 
dominant charge carries with S^, and contributing be¬ 
tween 50 % and 90 % of the total charge. While HCO^ is the 
dominant molecular ion, it contributes at most ~ 20 % of the 


® We define the molecular layer of a molecule such that the col¬ 
umn density of that molecule contained in it is equal to 90 % of that 
molecule’s total column density. It is centred at the position which has 
a largest fractional contribution to the total column density. 


charge in the inner regions and is severely depleted in the outer 
disk, R > 200 AU. Thus, while HCO^ being the most dominant 
molecular ion in the disk, even in the molecular layer its contri¬ 
bution to total charge is dwarfed by that of atomic ions such as 
C+, H+ and S+. 

This additional source of ionization not considered in the 
steady state approximation can contribute to the difference in 
observed values of x(e ). Furthermore, it is surprising that ob¬ 
servations sensitive to only a small region in the disk are able 
to recover disk average values relatively well. This is due to 
the density gradient towards the midplane; the disk average will 
draw heavily from values closer to the midplane. 


5. Summary 


In this paper Plateau de Bure Interferometer observations of the 
abundant molecules HCO+, J = (3-2), (1-0) and DCO+, J = (3- 
2) of DM Tau in parallel with a suit of thermo-chemical models 
provided the framework with which to study deuterium fraction¬ 
ation and the ionization fraction in DM Tau. 

Using combined -minimization and MCMC fitting tech¬ 
niques, we fitted a parametric model to the observations. HCO^ 
was found to exhibit a complex emission structure: (3-2) emis¬ 
sion had a peak intensity at r 50 AU and was considerably 
less extended than the (1-0) emission. HCO^ (1-0) emission was 
found to be co-spatial with CO emission ( Pietu et al.|2007| l. An 
inner hole of r « 50 AU in HCO^ is needed to recreate the obser¬ 
vations. Simultaneously fitting the (3-2) and (1-0) lines required 
the assumption that both lines had the same excitation temper¬ 
ature. DCO^ was also tentatively found to peak at r ^ 70 AU, 
consistent with the CO snowline, however higher resolution ob¬ 
servations are required to confirm this. 

We find Ru = A(HCO+)/A(DCO+) varies from 0.1 - 0.2 
between 50 and 430 AU, values considerably higher than both 
the cosmic abundance ~ 10“^ and those found in prestellar 
cores 0.035 ± 0.015. This is indicative of continued fractiona¬ 
tion throughout the disk lifetime. Both TW Hya and HD 163296 
exhibit similar high levels of deuteration, the later peaking at 
Ru ~ 0.3 ( |Qi et akpOOfil [Mathews et al.|2013| l. 
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Through chemical modeling we have also explored other 
physical parameters which can affect Rr,. The most influential 
parameters are the level of interstellar UV and the X-ray lumi¬ 
nosity of the central star, the dominant ionization sources in the 
molecular layer. X-rays interact with Rr, in a relatively straight 
forward manner: higher luminosities mean increased abundances 
of HCO^ and DCO^ and conversely for reduced luminosities. 
On the other hand, interstellar UV is more complex as it does 
not directly affect HCO^ and DCO^ but rather CO abundances. 
It enhances the DCO^ abundance across the entire disk, while 
suppresses the formation of HCO^ in the outer, r > 50 AU, disk. 
These effects can be disentangled through the dependence of Ro 
on radius with the later producing a more radially constant value 
of/Jo- 


Assuming a steady state system we estimate the electron 

10~^, consis- 


Dutrey et al. (2007 1 andjOberg et al. 


fraction of the HCO^ mole cular layer to be x(efl 
tent with lower limits from 
( |201 lb[ ). This values is high enough to induce MRI turbulence. 
An analysis of the dominant charge carries in the molecular layer 
show that HCO^ is the most dominant molecular ion, however 
atomic ions are considerably more dominant in all regions of 
the disk. Thus constraints on ionization from the abundance of 
HCO^ must take this into account. 
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